Dynamic charge correlations near the Peierls transition 
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The quantum phase transition between a repulsive Luttinger liquid and an insulating Peierls state 
is studied in the framework of the one-dimensional spinless Holstein model. We focus on the adiabatic 
regime but include the full quantum dynamics of the phonons. Using continuous-time quantum 
Monte Carlo simulations, we track in particular the dynamic charge structure factor and the single- 
particle spectrum across the transition. With increasing electron-phonon coupling, the dynamic 
charge structure factor reveals the emergence of a charge gap, and a clear signature of phonon 
softening at the zone boundary. The single-particle spectral function evolves continuously across 
the transition. Hybridization of the charge and phonon modes of the Luttinger liquid description 
leads to two modes, one of which corresponds to the coherent polaron band. This band acquires 
a gap upon entering the Peierls phase, whereas the other mode constitutes the incoherent, high- 
energy spectrum with backfolded shadow bands. Coherent polaronic motion is a direct consequence 
of quantum lattice fluctuations. In the strong-coupling regime, the spectrum is described by the 
static, mean- field limit. Importantly, whereas finite electron density in general leads to screening of 
polaron effects, the latter reappear at half filling due to charge ordering and lattice dimerization. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.30.+h, 71.45.Lr, 71.10.Pm, 02.70.Ss 



I. INTRODUCTION 

The Peierls transition is a hallmark feature of quasi- 
one-dimensional materials, such as conjugated polymers 
or organic charge-transfer salts. Peierls originally pre- 
dicted a band insulator for any finite electron-phonon 
coupling for a half-filled system due to an instability in 
the presence of perfect nestingii A full understanding of 
the physics requires to account for both quantum lat- 
tice fluctuations and strong electronic correlations)^ For 
example, it has been found theoretically that quantum 
fluctuations of the lattice give rise to a metallic state be- 
low a flnite critical coupling strength.— 

A minimal yet rich theoretical setting studied is the 
spinless Holstein model^ at half filling. With increas- 
ing electron-phonon coupling, a quantum phase transi- 
tion from a metal to a Peierls insulator (PI) with charge- 
density-wave order occursi^ The phase diagram has 
been mapped out with high a accuracyi^— The char- 
acter of the phase transition depends on the phonon 
frequencyi^ In the adiabatic regime considered here, the 
Peierls state is a band insulator, and the transition is 
accompanied by the softening of the phonon mode at 
the zone boundary. The critical coupling is rather small, 
so that our Monte Carlo method becomes very efficient 
(sec Sec. imp . In the nonadiabatic regime, small polarons 
form a polaronic superlattice, and the bare phonon mode 
hardens on approaching the transition)^ 

A metallic phase of one-dimensional spinless fermions 
is described by Luttinger liquid (LL) theory, with an in- 
teraction parameter Kp, even in the presence of a cou- 
pling to phonons )i2iii Until recently, from scaling analy- 
sis of the ground-state energy,—"— there was evidence for 
a crossover from an attractive LL {Kp > 1) at low phonon 
frequencies, to a repulsive LL {Kp < 1) at high phonon 



frequencies. However, large-scale density-matrix renor- 
malization group (DMRG) calculations of the charge 
structure factor have revealed that the LL is always re- 
pulsive, i.e. Kp < Ir^ The Peierls quantum phase tran- 
sition in the spinless Holstein model is of the Kosterlitz- 
Thouless (KT) type for all values of the phonon fre- 
quency, with Kp ~ 1/2 at the transition point )i^ The KT 
transition was previously known to occur in the nona- 
diabatic regime where the model can be mapped onto 
the XXZ Hamiltonian)^!^!^ The suppression of super- 
conducting correlations is ascribed to the spinless nature 
of the charge carriers, which prevents onsite pairing. In 
the spinful model at quarter filling, a transition from a 
LL to a (metallic) Luther-Emery liquid with finite spin 
gap but zero charge gap was discovcrcd)^^ An LL theory 
for gaplcss Peierls states has been dcveloped)^^ 

Given the continuous nature of the Peierls transition, 
one may argue that LL features should also survive at 
moderate couplings in the PI state, i.e. sufficiently close 
to the transition. Such a scenario will indeed emerge from 
our results for the single-particle spectral function, which 
reveal mixed charge and phonon modes (i.e., polarons) in 
both phases, with a finite single-particle gap in the Peierls 
phase. For the metallic phase, we shall compare to exact 
bosonization results for the spectral function)ii 

Dynamic properties of the spinless Holstein model 
have been studied numerically and analyticallyj^ii^^iSS 
(we refer to these papers for a more comprehensive re- 
view of work on the spinless Holstein model) except for 
the dynamic charge structure factor. The spinless Hol- 
stein model may be considered as the strong-Hubbard- 
U limit of the Hubbard-Holstein model, so that the 
formation of singlet bipolarons is suppressed. Away 
from half filling, it captures polaronic effects believed to 
be crucial in strongly correlated materials such as one- 
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dimensional MX chainsj^Zr— or quasi-two respectively 
three-dimensional manganiteSf^i^i^ 

Here we apply an exact continuous-time quantum 
Monte Carlo (CTQMC) method to revisit the Peierls 
transition in the spinless Holstein model. Apart from 
presenting a highly non-trivial test of the method for 
lattice problems, we make several contributions to the 
understanding of the physics near the transition. In par- 
ticular, we compute the dynamic charge structure factor, 
and relate the single-particle spectral function to LL the- 
ory and polaron physics on both sides of the transition, 
and to mean-field theory in the strong-coupling regime. 

The paper is organized as follows. We introduce the 
model in Sec. HIl and provide some details of the method 
for the present problem in Sec. IIIII A discussion of our 
results is given in Sec. |IVl and we conclude in Sec. El 

II. MODEL 

The Hamiltonian H = Hq + Hi takes the form 
^" = -*E(^Is +H-C-) -/^^n., (1) 

i ^ ^ i 

The first term Hq describes hopping of spinless fermions 
between neighboring sites with overlap integral t and free 
dispersion e{q) — ~2tcosq^ fi is the chemical potential, 
and the charge density operator is defined as hi = clc^ 
with eigenvalues 0, 1. The first part of Hi describes the 
free dynamics of the lattice in the harmonic approxima- 
tion; luq is the phonon frequency, and K = loqM. The 
second part couples the electron density to the lattice 
displacement. The half filled band (density n = 0.5) cor- 
responds to /i = 0, and we work in one dimension. 

Due to the use of first quantization (more convenient 
for the present method), the notation, in particular the 
coupling 7, differs from that used in previous work we 
shall compare to^i^ where the coupling term is written 

in the form —gijJoJ2ii^l + ~ i)- To reconcile the 

two notations, we use the dimensionless coupling con- 
stant A = 2ep/W, where Sp is the polaron binding en- 
ergy, and W — At IS the bandwidth; A is the relevant 
dimensionless ratio in the adiabatic regime wq < t, sepa- 
rating weak (A ^ 1) and strong coupling (A 1); small- 
polaron formation occurs at A ~ 1."^^ In terms of g we 
have £p = g^ujQ and A = g'^uJo/{2t), whereas in terms 
of 7 the relations are £p = '-f^/{2K) = j^/{2ujq) and 
A ~ '-f^/{ujQW). We use t as the unit of energy, and set 
h = kB = M = 1. 

III. METHOD 

The CTQMC is based on an exact diagrammatic 
expansion of the partition function around the non- 



interacting limitj^ the expansion converges for any fi- 
nite fermionic system at finite temperature. The fact 
that Wick's theorem holds for each configuration allows 
for a simple calculation of one and two-particle fermionic 
Green's functions. As discussed elsewhere^^i^ updates 
take the form of the addition and removal of single ver- 
tices, and optionally hipping Ising spins. The method is 
exact also for strong coupling, but becomes numerically 
less efficient due to large matrix sizes. The numerical ef- 
fort scales with the cube of the average expansion order; 
the latter depends linearly on the system size N, inverse 
temperature /3 = 1/T and (effective) coupling strength. 

The extension of the CTQMC method to electron- 
phonon problems has been discussed before Most im- 
portantly, as the algorithm is action based^^ the path 
integral over the phonon degrees of freedom can be per- 
formed exactly for a wide range of problems^^ (see also 
Ref . 113) ■ This leads to a retarded electron-electron inter- 
action of the form 

Si = - ff dTdT'Tln.ir) - i]Z?(r - r')[n.(r') - i] , 
Jo Jo j 

(2) 

where D{t) is the phonon propagator of the Holstein 
model [diagonal in real-space due to the onsite interac- 
tion in Eq. ([T])], and nj(T) is a Grassmann bilinear. The 
result is a purely fermionic algorithm, with no additional 
sampling of phonon degrees of freedom. The vertices 
have equal real space coordinates but in general different 
times t,t', distributed ^ D{t — t'). The range of the 
retarded interaction along the imaginary time is of the 
order but we find the algorithm to be efficient even 

for small loq. The weight to add a vertex is proportional 
to A/37V. In the present case, the spinless nature of the 
model makes only one spin direction appear in the algo- 
rithm. There is no sign problem, and we consider system 
sizes N = Al + 2 with periodic boundary conditions. 

Compared to previous QMC studies of the spinless Hol- 
stein model^i^i^i22i36 ^hc CTQMC algorithm is free of 
errors associated with a Trotter discretization. There is 
also no cutoff for the phonon Hilbert space, as the lattice 
degrees of freedom are integrated out exactly. 

Two key quantities for studying the Peierls transition 
arc the static charge structure factor 

Sp{q) - " ' (3) 

r 

and the dynamic charge structure factor given in the 
Lehmann representation as 

5p(g,^) = l^|(n|p,|m)|V^^™ 

mn 

X 5{E„, ~Er,~uj), (4) 

where \n) is an eigenstate with energy £"„, and we have 
defined = Y.r e'"'("r " {nr))/VN- 
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FIG. 1: (Color online) Density-density correlations Sp{r) = 
(hrho) (a) and charge structure factor (b). The coupling A = 
0.5 corresponds to the metallic Luttinger liquid (LL) phase, 
whereas A = 1.0 is in the insulating Peieris phase (PI), (c) 
Closeup of (b). Here /3t = N = 34, ujo = OAt and n = 0.5 
(half filling) . Lines are guides for the eye. 

We also consider the single-particle spectral function, 
Aiq,u^) = 1^1 {n\ c, \m) |^(e-/^^" + e"''^'") 

mn 

X 6{E^m -En-Lo). (5) 

For the continuation from imaginary time to real frequen- 
cies, we use the stochastic Maximum Entropy method.— 

IV. RESULTS 

We anchor our choice of parameters to previous 
work)2ii^ focusing on the adiabatic regime ujq < t (typi- 
cally realized in experiments) where the Peieris transition 
falls into the weak or intermediate coupling regime. In 
the latter, the CTQMC is by its nature particularly ef- 
ficient due to low average expansion orders Explicitly, 
following Ref. 0, we take ujq ~ OAt. We mainly con- 
sider two characteristic values of the coupling strength, 
namely A = 0.5 and 1 (equivalent to = 2.5 respectively 
5). The latter lie just above and below the Peieris transi- 
tion, see Fig. 1 in Ref. with the critical coupling known 
from DMRG calculations^ as Ac « 0.7 {g^ w 3_^. The 
same parameters have been considered in Refs. [20||23H23 . 
We shall also compare to the case ujq = O.lt considered 



in Ref. [l^, for which Ac « 0.4. The inverse temperature 
is set to I3t = N. 



A. Static charge correlations 

In Fig. [IJa) we show the density-density correlator 
Sp{r) ~ {firfiQ). The T = bosonization result isM 

K A 

(fixfio) = ~ r, 2^2 + ~^2ir cos{2kFx) , (6) 

with fcp ~ mr for the spinless case, equal to tt/2 for half 
filling. In the metallic phase, the LL parameter can be 
extracted from 

Kp^2Tr\imSp{q)/q. (7) 

For simplicity, we take Kp sa 2'KSp{qi)/qi, where qi is 
the smallest nonzero wavevector for a given lattice size. 

In the metallic LL phase, the amplitude of density cor- 
relations shows a power-law decay in real space. Since 
Kp = 0.96(1) < 1 for A = 0.5, we see a cusp in Sp{q) at 
q — 2kp. The linear form of Sp{q) at small q [Fig. [IJc)] 
is evidence for a metallic state at A = 0.5. 

In the Peieris phase (A = 1), the density correlation 
functions display long-range order (a charge density wave 
with one fermion on every other site, and a staggered 
lattice displacement with periodicity tt). The small-g be- 
havior of Sp{q) becomes nonlinear, see Fig.[TJc). Via the 
continuity equation, the long- wavelength limit of Sp{q) 
can be directly related to the Drude part of the opti- 
cal conductivityfi^ and the absence of linear behavior of 
Sp{q) as g — > is equivalent to a vanishing Drude weight 
in the insulating Peieris state. As a result of long-range 
charge-density- wave order, Sp{q) diverges with system 
size at 5 = TT = 2fcF, and serves as an order parame- 
ter. The Peieris state is quickly suppressed away from 
the commensurate density n = 0.5i^ 

The LL parameter Kp has been calculated as a func- 
tion of g with high accuracy in a large-scale DMRG study 
including finite-size scaling with up to 256 sites.— As the 
CTQMC is not capable of studying such large systems, 
and due to the additional complication of finite tempera- 
tures in our method [Eqs. ^ and ([7]) only hold at T = 0] 
we do not attempt a quantitative comparison here. In- 
stead, we have checked for both ujq — 1 and 0.1 that 
even on small systems {N ~ 22) the CTQMC correctly 
reproduces Kp < 1 (repulsive LL), albeit with signifi- 
cant finite-size effects near gc as expected for a KT phase 
transitions^ 



B. Excitation spectra in the metalHc phase 

We now turn to the discussion of the dynamical cor- 
relation functions, namely the dynamic charge structure 
factor and the single-particle spectral function. To our 
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FIG. 2: (Color online) (a) Dynamic charge structure factor 
Sp{q,uj) and (b) single-particle spectral function A{q,Ltj) in 
the metallic LL phase. Horizontal lines indicate (a) the bare 
phonon frequency ojo and (b) the chemical potential fi = 
(long dash) and ±oJo (short dash). We have limited the inten- 
sity range to reveal the subdominant features; the full range 
for (a) is shown in Fig. |6jb); for (b) see also Fig. [S] Parame- 
ters are ujq = OAt, X = 0.5, and pt = N = 34. 



knowledge, the former has not been calculated for the 
present model before. In all intensity plots, we have re- 
frained from using any kind of interpolation. This makes 
the limited (yet significantly higher than in previous 
numerical work) momentum resolution well discernible. 
The energy resolution was taken to be O.Olt or better. 

The spectra in the LL phase (A = 0.5) are shown in 
Fig. [5] The dynamic charge structure factor Sp{q,uj) 
[Fig. [2{a)] shows the familiar continuum of particle-hole 
excitations of fermions in one dimension. At long wave- 
lengths, wc have a linear, gapless mode with disper- 
sion N{q,uj) w N{q)5{uj — Vpq) expected for the metallic 
phase. Near the zone boundary, the width of the con- 
tinuum approaches the bare bandwidth W ^ At. In 
addition, because of the coupling between the electron 
density and the phonons in the Hamiltonian ([1]), a sig- 
nature of the renormalized optical phonon mode is vis- 
ible at low energies. Already below the Peierls transi- 
tion at Ac ~ 0.7, it shows partial softening near q = n. 
Our findings are in excellent agreement with previous 
calculations of the phonon spectral functiouj^i^ see also 



FIG. 3: (Color online) Single-particle spectral function 
A{q,ijj) in the metallic LL phase. Parameters as in Fig. [2jb). 



Figs, (nja) and (b) discussed below. Hence, despite the 
nondispersive nature of the Einstein phonons in Eq. ([1]) 
the coupling gives rise to a momentum dependence of the 
renormalized phonon frequency. 

The single-particle spectrum A{q, lo) in the LL phase 
[Fig. [DJb)] reveals a dominant band (note the loga- 
rithmic scale in all our intensity plots) closely track- 
ing the cosine dispersion of the noninteracting problem, 
e(g). Close to the Fermi level, inside the coherent in- 
terval id e [— woi'^o]- the main peak becomes noticeably 
sharper 1^ As noted previously^^'^ the incoherent part 
(1^1 > Wo) of the spectrum is broadened due to multi- 
phonon processes;^ the width is roughly = 2.5. 

As expected, but not discussed in the numerical 
literature^ii^ (where the spectra are mostly interpreted 
in terms of polaron theory), the spectrum shares the 
main features of the exact bosonization result for a lin- 
ear fermion dispersionJ^ For comparison, we show the 
results in Fig.[5Jb) again in Fig.[3]using a representation 
similar to that in Ref. [ill The bare charge and phonon 
modes of the LL Hamiltonian of the spinless model hy- 
bridize and broaden; see also Ref. [l6l This is equivalent 
to the formation of coherent polaronic quasiparticles, i.e. 
electrons dressed with phonons, with increased effective 
mass. At fcp = tt/2 (we focus on cj < 0), there is a domi- 
nant zero-energy peak and a phonon satellite with thresh- 
old energy —loq whose width depends on the coupling 
strength. With q moving away from fcp, spectral weight 
is transferred to the "phonon peak" (in the nomencla- 
ture of Meden et alJ^) which eventually becomes the 
dominant excitation and follows e{q). Since in contrast 
to Ref. flU the bandwidth W ^ coq, we observe multiple 
phonon satellites which merge into a broad band. This 
picture is consistent with previous numerical^ii^ and an- 
alytical results^i (though difficult to see in the latter) in 
the metallic phase. Similar results and an explicit com- 
parison to LL theory have been reported for the spinful 



5 



0.01 



0.1 



10 



0.01 



0.1 



10 



(a) Sp(q,co) 








F 






(c) A(q,(o) A 



71/2 

q 



FIG. 4: (Color online) (a) Dynamic charge structure factor 
Sp{q, uj) and (b) single-particle spectral function A{q, lj) in the 
Peierls phase, (c) Closeup of (b). Horizontal lines indicate (a) 
the bare phonon frequency tJo, (b) the chemical potential ^ = 
(long dashed) and iwo (short dashed), and (c) ±ajo. The 
solid line in (c) is the shifted single-polaron band dispersion 
ep(g) + A/2 for the same parameters,— mapped from [0, vr] to 
[fcp , tt] to facilitate comparison with half filling. The double- 
headed arrows in (a) and (b) designate the mean-field gap 
Amf = Ep = 2t. The single-headed arrows in (a) and (c) 
indicate the Peierls gap A « 0.75£. In (a) we have limited 
the intensity range to reveal the subdominant features; the 
full range is shown in Fig. Oc). Parameters are ojo = QAt, 
A = 1.0, and pt ^ N ^ 34. 



C. Excitation spectra in the Peierls phase 

In Fig. 0] we report the dynamic charge structure factor 
and the single-particle spectrum in the Peierls state, at 
A = 1. Mean-field results for A{q,u}) to be compared to 
Fig. mb) are shown in Fig. [SJ 

In the case of Sp{q,uj), see Fig. HJ^a), we note that al- 
most the entire spectral weight is located close to a; 0, 
g = TT. All other features, including the particle-hole con- 
tinuum, have a spectral weight that is at least an order 




FIG. 5: (Color online) Mean-field result for the single-particle 
spectrum A{q,Lj) of the adiabatic (ij;o — 0), spinless Holstein 
model at T = 0. Here A^f = 2t (arrow); we have replaced 
the S peaks by Lorentzians with scaling parameter 0.3t. 



of magnitude lower (see caption). The accumulation of 
spectral weight at uj — 0, q ^ tt ^ 2kF is a result of the 
complete softening of the renormalized phonon dispersion 
and is characteristic of the displacive Peierls transitioni^ 

The renormalization of the phonon mode with increas- 
ing A is well visible in Fig. |6l which shows the low- 
energy range of Sp{q,uj) on a larger scale. Already for 
A ~ 0.25 [Fig. inja)] a partial softening is visible close to 
q = TT. This trend continues with increasing A, until the 
phonon has become gapless at g = tt in the Peierls state 
[Fig. [SJc)]. The softening scenario agrees well with pre- 
vious worki^ i^^i^^i^^ From Sp{q,uj) we cannot detect the 
previously reported splitting of the phonon mode into two 
branches, one of which hardens toward luq with increas- 
ing X-^^ this may be a result of the very small spectral 
weight of the upper branch in the adiabatic regime. 

The linear long-wavelength mode characteristic of the 
metallic phase [Fig. dja)] is completely suppressed with 
increasing A, see Fig. |6l From Fig. IH^a), we can in- 
fer a gap of about 0.75i (left arrow) above which the 
particle-hole continuum starts. However, we can also es- 
timate a second characteristic energy of about 2t (right 
arrow) where the spectral weight of the continuum in- 
creases noticeably. These two energy scales are related 
to the single-particle spectrum below. Our results for the 
dynamic charge structure factor are in accordance with 
previous calculations of the phonon spectral function for 
the present modeli^i^i By comparison, the Luther-Emery 
phase of the spinful Holstein model is characterized by a 
soft phonon mode in addition to the particle-hole contin- 
uum which remains gapless for g — > OJ^ 

The single-particle spectral function in the Peierls state 
is shown in Fig. IDJb). As hinted at by the results for 
Sp{q,uj), it is found to consist of two sets of features. 
There is a main, high-energy band that follows the free 
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FIG. 6: (Color online) Dynamic charge structure factor 
Sp{q,uj) at low energies. Horizontal lines indicate the bare 
phonon frequency loq. Here ljo = OAt, pt — N — 22 and (a) 
A = 0.25 (LL), (b) A = 0.5 (LL), and (c) A = 1 (PI). 



dispersion far away from fcp, and is split by the mean- 
field gap (see below) Amf ^ Sp — 2t (indicated by the ar- 
row). This value agrees with the gap 2t discussed above. 
The main band also reveals backfolded shadow bands for 
q > kp and w < wq respectively for q < kp and uj > uiq. 
Additionally, within the mean-field gap, we have disper- 
sive low-energy modes with a smaller gap A « 0.75t at 
kp. We refer to this actual gap (determining the low- 
energy properties of the system) as the Peierls gap. 

To first understand the high-energy features of A{q, lu), 
we consider the adiabatic limit wq = of the spinless 
Holstein model,— 



H' = Ho- i/v ^ - 7 E - • 



(8) 



With the mean-field ansatz Qi = (—1)'^, we obtain the 
bands E±{q) = ±[e{q)'^ + A^f/4]i/2 ^-^^^ g, gap A^f at 
kp = 7r/2. The mean-field gap is related to the lattice 
order parameter, Amf ~ and as for displaced oscilla- 
tor states the lattice shift on occupied sites is <^ ~ y/e^ so 
that we expect Amt ~ £p, m agreement with the QMC 
results. Rather than solving the gap equation, we simply 
set Anif equal to the gap in our numerical data. 

The resulting single-particle spectral function^ con- 
tains two branches following E±[q), centered around 
q = Q respectively tt, with spectral weights 'w±(q) = 
[1 ± e{q)/\E±{q)\]/2, and is shown in Fig. [3 Setting 
Amf = Ep and replacing the 5 functions by Lorentzians 
with scaling parameter 0.3i, we find excellent agree- 
ment with the high-energy features in Fig. |31[b) . In par- 
ticular, the backfolded shadow bands (including the q- 




FIG. 7: (Color online) (a) Dynamic charge structure factor 
Sp{q, uj) and (b) single-particle spectral function A[q, cj) in the 
Peierls phase. Horizontal lines indicate (a) the bare phonon 
frequency wo, (b) the chemical potential ^ = 0. Arrows des- 
ignate the mean-field gap Amf = Cp = l.&t. In (a), we have 
limited the intensity range to reveal the subdominant features. 
Parameters are ojo = 0.14, A = 0.8, and pt = N = 22. 



dependence of the spectral weight) emerge as a natu- 
ral feature of gapped systems with competing periodic 
potentials4^ The fact that the high-energy features of 
the quantum case ujq > are well captured by mean- 
field theory justifies the notation A,„f introduced above. 

Let us now focus on the low-energy features inside the 
mean- field gap [Fig. IDJc)]. The absence of the latter in 
the mean-field results for the static limit suggests that 
they are related to quantum phonon fluctuations. In- 
deed, a particle added to the half-filled Peierls state [as 
described by A{q > kp,uj > 0), cf. Eq. ([5])] is enabled 
by lattice fluctuations to propagate as a polaron. For 
the present parameters, the dimerization of the lattice 
opens a small but flnitc gap A ^ A„if in the correspond- 
ing polaronic band. Support for this picture comes from 
several directions. First, the width of the polaron bands 
in Fig. |3^c) is about 0.2t, in excellent agreement with 
the single-polaron results for the same parameters.— It 
also shows the well known flattening of the dispersion 
with increasing gj^ii^ To demonstrate this agreement, 
we include in Fig. Ul^c) the band dispersion ep(g) of a 
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single polaron for the same parameters 4i To compare 
to half filling, we shift the dispersion by A/2, and map 
it from [0, tt] to [feFiTr]. This yields excellent agreement 
for all q. Second, these bands have large electronic spec- 
tral weight near q = kp (corresponding to q = for 
a single polaron) but very small weight far away from 
fcp; the character of the coherent quasi-particle changes 
from electronic to phononic when the band energy in- 
tersects the bare phonon frequency ; ^^i*^^ see also Fig. |3l 
Third, the energy gap A of the polaron bands, which 
is a property of the many-particle Peierls state, matches 
the corresponding excitation in the phonon spectral func- 
tion for the same parameters?^ and agrees well with the 
phonon signature in Fig. |6l^c) near q = 7r/2. For a sin- 
gle polaron, no phonon renormalization occurs, and the 
polaron band dispersion is clearly visible in the phonon 
spectral function.— The spectral function for the same 
parameters as Fig. 2] has been calculated by cluster per- 
turbation theoryj^ Whereas the polaron bands are diffi- 
cult to identify in these approximate results, their sig- 
nature inside the mean-field gap is clearly visible in the 
exact density of states<^ The QMC results are in ex- 
cellent agreement with the T = wavevector-resolved 
spectral function for these parameters obtained by ex- 
act diagonalization4i An alternative explanation for the 
polaron bands in terms of thermal excitations based on 
QMC simulations at much higher temperature was given 
beforci^ Polaron bands inside the static mean-field gap 
were also predicted analytically?^ and may be related to 
soliton excitations*^ 

From LL theory, the existence of the low-energy po- 
laron modes is a necessary consequence of the continuous 
nature of the Peierls transition. The hybridized charge 
and phonon modes of the metallic phase evolve continu- 
ously with the coupling strength A. One of the resulting 
two modes represents the backfolded high-energy band, 
whereas the other mode is the polaron band. This struc- 
ture is already visible in the LL phase, see Figs, ^h) 
and El The Peierls gap A in the polaron band opens 
exponentially slowly at the KT transition. 

In the light of previous work on many-polaron systems, 
the validity of the single-polaron picture at half filling 
is not obvious. Away from half filling, the spinless Hol- 
stein model is metallic, and the overlap of the (extended) 
lattice distortions of individual carriers has been found 
to lead to a renormalization toward weak coupling and 
a loss of polaronic signatures j^i^i^ However, the elec- 
trons constituting the insulating Peierls state are ordered 
in the dimerized lattice potential. Hence, an additional 
particle or hole does not undergo screening and behaves 
as a polaron, which moves in a dimerized but fluctuat- 
ing potential. The most obvious impact of this Peierls 
background is the opening of the gap A. 

To complete the understanding of the dynamic charge 
correlations, we consider the strong-coupling regime. 
From single-polaron theory, we expect a strong reduc- 
tion of the polaron bandwidth and the electronic spec- 
tral weight; the latter also becomes almost independent 



of qJ^t^ The many-body Peierls gap A should increase 
with increasing dimerization (~ VA) of the lattice com- 
pared to phonon fluctuations (~ -^/wo), leading to a sup- 
pression of the coherent polaron motion. 

Exact T = results for the single-particle spectrum in 
the Peierls phase have been reported by Sykora et al.^ 
for N = 8, ujQ = O.lt and A = 0.8. On the linear scale 
of their Fig. 4, no features of the kind discussed here are 
visible. These results fall into the strong-coupling regime 
with respect to the Peierls transition, A ~ 2 Ac. In par- 
ticular, the mean-field gap A^f 3> wq- As demonstrated 
by the QMC results for the parameters of Sykora et al 
shown in Fig. [71[b), the spectral weight of the polaron 
band is extremely small. At fcp, the gap A almost co- 
incides with the mean-field gap Amf- The spectrum is 
dominated by the mean-field features of Fig. [51 As a 
result, the dynamic charge structure factor in Fig. [7Ja) 
reveals only one energy scale (A^f w A) for the onset of 
the particle-hole continuum, in addition to the signatures 
of phonon softening at low energies. A similar result is 
expected for = OAt at much stronger coupling where 
the CTQMC method becomes inefficient. 

Finally, the features of the spectral function in the 
Peierls state discussed here are specific to the adiabatic 
regime < t. In the nonadiabatic case wq ^ t, the crit- 
ical coupling for the Peierls transition is large, so that 
electrons become heavy, small polarons already below 
Ac ^ The corresponding polaron bands are extremely nar- 
row and have very small electronic spectral weight J^i^ 
Similarly, the dispersion of the high-energy features is 
also substantially reduced. 



V. CONCLUSIONS 

We have applied the exact continuous-time Monte 
Carlo method to the spinless Holstein model in the adi- 
abatic regime. We find that the method is well suited to 
study the Peierls metal-insulator transition on the lattice. 
We have complemented previous work by computing the 
dynamic charge structure factor, and showing that is 
allows to track the softening of the phonon mode and 
the appearance of a charge gap with increasing electron- 
phonon coupling. For the single-particle spectral func- 
tion, we explicitly demonstrated that the mixed charge 
and phonon modes of Luttinger liquid theory constitute 
the spectrum both below and above the Peierls tran- 
sition. Most importantly, the spectrum in the Peierls 
phase consists in general of high-energy features readily 
understood in the static limit, and gapped low-energy 
polaron bands originating from quantum lattice fluctu- 
ations. These modes disappear in the strong-coupling 
regime. Our findings reconcile and unify previous and 
present results with both polaron theory and Luttinger 
liquid theory. 

The spinless Holstein model captures two essential 
consequences of electron-phonon coupling, namely the 
formation of polarons (dressed quasi-particles) and the 
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Peierls metal-insulator transition at commensurate fill- 
ing. For any finite electron-phonon coupling, the elec- 
trons acquire some polaronic character. Based on this 
and previous work, we can conclude that starting from 
the low-density limit the well-defined polaron signa- 
tures known from the single-electron case^^ first be- 
come washed out with increasing density in the metallic 
phascj^ However, once we approach the Peierls phase at 
half filling, mutual screening is strongly suppressed due to 
charge order and lattice dimerization, and clear polaron 



signatures reappear in the single-particle spectrum. 
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